tr"  53  CO 


AD-A149  313 


TECHNICAL  REPORT  BRL-TR-2598 


AVERAGING  EFFECTS  IN  MODELS  OF 
THREE-DIMENSIONAL  TWO-PHASE  FLOWS 

Aivars  K.R.  Celmins 

October  1984 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


US  ARMY  BALLISTIC  RESEARCH  LABORATORY 

ABERDEEN  PROVING  GROUND.  MARYLAND 


11  O7 


00  S 


84 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  it  to  the  originator. 


Additional  copies  of  this  report  may  be  obtained 
from  the  National  Technical  Information  Service, 

U.  S.  Department  of  Commerce,  Springfield,  Virginia 


22161. 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position,  unless  so  designated  by  other 
authorized  documents. 


The  use  of  trade  names  or  manufacturers’  names  in  this  report 
does  not  constitute  indorsement  of  any  commercial  product. 


;LASS:FICATlON  OF  This  PAGE  (When  Dele  Entered) 


REPORT  DOCUMENTATION  PAGE 


4.  TITLE  fond  SubtlU«) 

Averaging  Effects  in  Models  of  Three-Dimensional 
Two-Phase  Flows 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  recipient's  catalog  number 


5.  TYPE  OF  REPORT  &  PERIOD  COVERED 


6.  PERFORMING  ORG.  REPORT  NUMBER 


AUTHORS) 

Aivars  K.  R.  Celmins 


B.  CONTRACT  OR  GRANT  NUMBERfa) 


9  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

PS  Army  Ballistic  Research  Laboratory 

ATTN:  AMXBR-IBD 

Aberdeen  Proving  Ground,  Maryland  21005-5066 


10.  PROGRAM  ELEMENT.  PROJECT.  T  ASK 
AREA  &  WORK  UNIT  NUMBERS 


12.  REPORT  DATE 


11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

US  Army  Ballistic  Research  Laboratory 
ATTN:  AMXBR-OD-ST 

Aberdeen  Proving  Ground,  Maryland  21005-5066 


4  MONITORING  AGENCY  NAME  »  AOORESSf//  dllterent  Irom  Controlling  Olllce)  15.  SECURITY  CLASS,  (o I  thle  raporf) 


October  1984 


13.  NUMBER  OF  PAGES 

82 


16  DISTRIBUTION  STATEMENT  (o(  thle  Report) 


UNCLASSIFIED 


1S«.  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 


Approved  for  public  release;  distribution  is  unlimited. 


17.  DISTRIBUTION  STATEMENT  (at  the  abetted  entered  In  Block  30,  II  dllterent  Iroai  Report) 


' 9  KEY  WORDS  (Continue  on  reveree  ride  It  neceeeery  end  Identity  by  block  number) 

Volume  Average  Interior  Ballistics  Calculations  Lattice  Arrangements 
Gas-Particle  Mixture  Undulations  of  Flow-Undulations  Mean  Particle  Distance 
Two-Phase  Flow  Particle  Induced  Undulations  Neighbor  Distance 
Gas  Volume  Fraction  Bounds  of  Undulations  Smoothing  by  Averaging 
Particle  Aggregate  Boundary 


20  ABSTRACT  (Continue  an  reveree  ride  II  rmceeeery  end  Identify  by  block  number) 

In  order  to  avoid  the  treatment  of  individual  particle  motion  in  two- 
phase  fLow  description,  one  can  use  volume  averaged  descriptors.  However, 
if  the  averaging  volume  is  too  small  then  individual  particles  can  cause 
large  undulations  of  the  averages.  In  this  report  the  amplitude  of  such 
undulations  is  estimated  by  a  bound.  It  is  shown  that  in  general,  one  must 
average  over  30-150  particles  in  order  to  obtain  reasonably  smooth  averages. 
A  consequence  for  interior  ballistics  calculations  is  that  volume  averages 


FORM 
I  JAM  73 


EDITION  OF  I  NOV  6S  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fNhwn  Dele  Entei-d) 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAOE(1Wi«n  Data  KntamQ 


20.  ABSTRACT  (Continued) 

may  be  used  to  model  the  core  flow  but  in  general  they  cannot  be  used  to 
resolve  radial  structures  In  two  phases.  Also  presented  are  examples  of  the 
representation  of  a  gas-particle  aggregate  boundary  In  terms  of  the  gas 
volume  fraction,  and  an  approximate  expression  is  derived  for  the  transition 
curve. 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGEfltTian  Data  Entarad) 


TABLE  OF  CONTENTS . 

LIST  OF  FIGURES . 

INTRODUCTION . 

PRINCIPAL  RESULTS . 

GAS  VOLUME  FRACTION  AS  A  FUNCTION  OF  THE  AVERAGING  VOLUK . . 

GAS  VOLUME  FRACTION  PROFILES . 

DISCUSSION  OF  THE  RESULTS . 

APPENDIX  A.  LATTICES . 

APPENDIX  B.  NEIGHBOR  DISTANCE  AND  JEAN  DISTANCE . 

APPENDIX  C.  ALGORITHM  FOR  GAS  VOLUME  FRACTION  CALCULATION 

APPENDIX  D.  GAS  VOLUJE  FRACTION  AT  A  PARTICLE  AGGREGATE 

BOUNDARY . 

APPENDIX  E.  UNDULATIONS  OF  AVERAGE  FUNCTIONS . 


DISTRIBUTION  LIST 


'.I  ■  !■» 


■  l  ■  \  1  L  ■  p.  ■  '■ 


.1  .'■'.'•IVL'l.A'  TTmr 


T 


LIST  OF  FIGURES 


’  Agure  Page 

1  Averaging  Sphere  Radius  for  Given  Tolerance  |^a|t0^ . 10 

2  Minimum  Number  of  Particles  in  an  Averaging  Volume 

for  given  Tolerance  | .^oi  1 1 o  1  an<^  ^  =  ^ 

*  3  Gas  Volume  Fraction  Dependence  on  Averaging  Radius. 

Square  Cylinder  Lattice,  a  =  0.5.. . 16 

i  4  Gas  Volume  Fraction  Dependence  on  Averaging  Radius. 

Squire  Cylinder  Lattice,  i  =  O.9.... . 17 

5  Gas  Volume  Fraction  Dependence  on  Averaging  Radius. 

Triangular  Cylinder  lattice,  i=  0.9....... . 18 

6  Gas  Volume  Fraction  Dependence  on  Averaging  Radius. 

Leap-Frog  Square  lattice,  ci=  0.9 . 19 

7  Gas  Volume  Fraction  Dependence_on  Averaging  Radius. 

Leap-Frog  Triangular  Lattice,  n=  0.9..... . . . 20 

3  Extremes  of  Gas  Volume  Fraction  as_Functlons  of  Averaging 

Radius.  Square  Cylinder  lattice,  o=  0.5 . 21 

9  Extremes  of  Gas  Volume  Fraction  as  Functions  of  Averaging 

Radius.  Triangular  Cylinder  Lattice,  i*  0.5 . 22 

10  Extremes  of  Gas  Volume  Fraction  as  Functions  of  Averaging 

Radius.  Leap-Frog  Square  lattice,  u=  0.5 . 23 

11  Extremes  of  Gas  Volume  Fraction  as  Functions  of  Averaging 

Radius.  leap-Frog  Triangular  lattice,  •*»  0.5 . 24 

12  Extremes  of  Gas  Volume  Fraction  as_Functions  of  Averaging 

Radius.  Square  Cylinder  Lattice,  i*  0.9. . 25 

13  Extremes  of  Gas  Volume  Fraction  as  Functions  of  Averaging 

Radius.  Triangular  Cylinder  Lattice,  <i=  0.9 . 26 

14  Extremes  of  Gas  Volume  Fraction  as  Functions  of  Averaging 

Radius.  leap-Frog  .Square  Lattice,  a*  0.9 . 27 

15  Extremes  of  Gas  Volume  Fraction  as  Functions  of  Averaging 

Radius.  Leap-Frog  Triangular  Lattice,  0.9.. . 28 

k 

16  Extreme  Deviations  of  Gas  Values  Fraction. 

leap-Frog  Triangular  Lattice,  a*  0.9 . . . 29 

17  Extreme  Deviations  of  Gas  Volume  Fraction. 

Leap-Frog  Triangular  Lattice,  «■  0.9 . 30 


_ • 


5 


LIST  OF  FIGURES  (Continued) 

Figures  Page 

18  Extreme  Deviations  of  Gas  Volume  Fraction  for  All 

Lattices.  a  m  0.5......... . .........31 

19  Extreme  Deviation  of  Gas  Volume  Fraction  for  All 

lattices.  a  =  0.9..................... . . . 32 

20  Gas  Volume  Fraction  at  a  Particle  Aggregate 

Boundary,  a  =  0.5..... . . . . . ..35 

21  Gas  Volume_Fraction  at  a  Particle  Aggregate 

Boundary,  a  3  0.7 . . . 37 

22  Gas  VolumeFraction  at  a  Particle  Aggregate 

Boundary,  u  3  0.9 . . . 38 

23  Extreme  Values  and  Estimated  Bounds  of  Gas  Volume  Fraction  at  a 

Particle  Aggregate  Boundary,  a  =  0.5,  R/I^  3  1 . 39 

24  Extreme  Values  and  Estimated  Bounds  of  Gas  Volume  Fraction 

at  a  Particle  Aggregate  Boundary,  a  3  0.9,  R/ 3  1 . 40 

25  Extreme  Values  and  Estimated  Bounds  of  Gas  Volume  Fraction 

at  a  Particle  Aggregate  Boundary.  ci=  0.9,  R/ L^  -  2 . 41 

A.l.  Square  Cylinder  lattice . 48 

A. 2.  Triangular  Cylinder  lattice.. . 48 

A. 3.  Leap-Frog  Square  Lattice........ . 50 

A. 4.  Leap-Frog  Triangular  lattice. . 50 

C.l.  Intersection  of  Spheres . . . 63 


1.  INTRODUCTION 


A  common  method  for  the  derivation  of  a  manageable  flow  description  and 
of  governing  equations  for  flows  of  gas-particle  mixtures  is  the  averaging 
of  local  flow  properties  over  a  volume.*  The  averaging  produces  from  the 
heterogeneous  local  flow  properties  smooth  flow  parameter  functions  which 
provide  representative  descriptions  of  average  properties  of  the  particle 
aggregate  and  of  the  gas  between  particles.  The  smoothing  is  most  effective 
if  the  averaging  volume  is  so  large  that  the  contributions  of  single 
particles  to  the  average  is  negligible.  Therefore,  it  is  reasonable  to 
choose  a  large  averaging  volume.  On  the  other  hand,  any  volume  averaging 
smoothes,  reduces  and  distorts  flow  structures  particularly  those  which  have 
an  extension  smaller  than  the  averaging  volume.  Hence,  in  order  not  to  lose 
flow  structures  of  interest  one  should  choose  a  small  averaging  volume.  In 
order  to  make,  under  these  conditions,  a  rational  choice  of  the  size  of  the 
averaging  volume  one  needs  a  quantitation  of  the  smoothing  effect  of  volume 
averaging. 

In  this  report  a  quantitation  of  the  smoothing  effect  is  obtained  from 
an  investigation  how  undulations  of  the  gas  volume  fraction  function 
oc  depends  on  the  size  and  distribution  of  particles  and  on  the  size  of  the 
averaging  volume.  Other  flow  parameters  can  be  shown  to  have  undulations 
that  are  proportional  to  the  undulations  of  a.  The  result  of  the 
investigation  is  an  estimate  of  bounds  for  the  undulations  of  a  in  terms  of 
the  averaging  parameters.  Using  this  estimate  one  can  choose  a  tolerance 
level  for  the  undulations  and  obtain  a  corresponding  minimum  size  of  the 
averaging  volume.  Flow  structures  with  extensions  smaller  than  the  chosen 
averaging  volume  are  distorted  and  reduced  by  the  averaging,  that  is,  they 
are  not  correctly  represented  by  the  average  functions. 

We  illustrate  the  smoothing  of  flow  structures  by  considering  a  plane 
boundary  of  a  region  with  uniform  particle  distribution.  In  concept,  such  a 
boundary  is  a  narrow  transition  zone  between  the  regions  with  a  =  1  (gas 
only)  and  a  =  a  (the  average  gas  volume  fraction  in  the  mixture  region). 
The  a  obtained  by  volume  averaging  has  Instead  a  relatively  wide  transition 
zone  which  is  spread  out  over  a  diameter  length  of  the  averaging  volume. 
The  width  of  the  transition  zone  cannot  be  reduced  arbitrarily  without 
penalty,  that  is,  without  increasing  undulations  of  the  average  flow 
parameters.  Hence  one  has  to  choose  between  a  smooth  representation  of  the 
average  flow  field  and  an  accurate  representation  of  the  boundary  of  the 
region.  The  estimated  bounds  of  the  undulations  can  help  one  to  make  the 
choice  rationally. 


*Ce£mim,  A ivaAA  K.R.  and  Schm^otX,  Jamu  A.,  ThAte-dweMionaJL  modeling  o ^ 
gcu>  -cumbtuting  ioEtd  two-pha, 6e  filom,  PA.oceeding-5  ofi  the  Tkoid  Uulti-phaAe 
FCuw  and  Heat  TfiaM^eA.  SympoAtum-Wosikihop,  pp.  681-698,  Ed.,  T.N.  Vzzvioglu, 
18-20  A pK  1983,  M iami  Beach,  FL. 
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In  Section  2  we  present  the  principal  results  of  the  investigation:  a 
quantitative  estimate  of  bounds  for  particle  induced  undulations  of  average 
flow  parameters.  Particulars  about  the  derivation  of  the  estimate  are  given 
in  Section  3 .  Section  4  provides  numerical  examples  of  the  representation 
of  a  particle  aggregate  boundary  in  terms  of  a  .  The  transition  profile  can 
be  computed  by  an  approximate  formula  which  is  also  given  in  Section  4  .  All 
calculations  are  for  a  spherical  averaging  volume  and  spherical  particles 
arranged  in  regular  three-dimensional  lattices.  However,  by  considering 
different  types  of  lattices  we  have  obtained  results  that  are  representative 
for  all  regular  particle  arrangements.  Also,  the  restriction  to  spherical 
particles  is  of  little  consequence  if  the  particles  are  small  compared  to 
the  averaging  volume.  Section  5  contains  a  discussion  of  the  results. 

2.  PRINCIPAL  RESULTS 

In  Section  3  we  consider  L^e  gas  volume  fraction  a  within  an  averaging 
sphere  with  the  radius  R,  and  investigate  undulations  of  a  due  to  the 
location  of  the  sphere.  Let  the  particles  be  arranged  in  a  regular  three- 
dimensional  array,  let  the  particle  radius  be  s,  and  let  a  be  the  limit  of 
the  gas  volume  fraction  as  the  averaging  sphere  becomes  infinite.  We  define 
a  mean  distance  1^  between  the  particle  centers  by 

Lju  =  2s(l  *  o)  _1/3  .  (2  .1) 

A  motivation  for  this  definition  is  given  in  Appendix  B,  where  it  is  also 
shown  that  Is  10-24%  larger  than  the  smallest  distance  between  particle 
centers .  If  the  averaging  sphere  is  finite,  then  one  obtains  Instead  of  a  a 
value  of  the  gas  volume  fraction  that  depends  on  the  position  of  the 
sphere.  Let  Ao  be  the  difference  between  an  actual  value  of  a  and  the  limit 
value  a : 


Aa  =  a  -  a  .  (2 .2) 


Aa  generally  depends  on  the  position  of  the  sphere  as  well  as  on  a  and 
L^/R.  In  Section  3  we  show  that  its  magnitude  is  bounded  by 

|Aa|  <  0.5  (1  -  a)  a2  (I^/R)2  (2.3) 

independently  of  the  position  of  the  averaging  sphere.  Eq .  (2.3)  is  based 
on  sample  calculations  with  different  lattices  for  1  <  R/L^  <  4  and  0.5 
<  at  <  0.9.  Because  the  formula  also  produces  the  correct  limit  Aa  «  0  for 
a  =  1,  it  can  be  used  as  an  estimate  for  all  a  >  0.5.  Extrapolation  to 


smaller  values  (a  <  0.5)  should  be  done  with  reservations,  and  the  same 
applies  to  extrapolations  to  R/  >  4.  The  domain  R/ <  1  is  of  little 
practical  interest  because  there  the  undulations  are  too  large. 

Eq.  (2.3)  is  graphically  displayed  in  Figure  1  as  a  relation  between  a 
maximum  tolerance  |^a|t0]^  and  the  corresponding  R/  which  guarantees  that 
the  undulations  are  less  than  the  tolerance.  The  three  curves  are  shown  for 
a  =  0.5,  0.67  and  0.9,  respectively.  The  value  a  =  0.67  is  according  to  Eq . 
(2.3)  the  worst  case,  that  is,  in  this  case  one  needs  the  largest  R/  I1n  in 
order  to  suppress  the  undulations.  The  solid  lines  indicate  the  domain  for 

which  the  formula  (2.5)  was  developed  and  tested  by  sample  calculations. 
Extrapolations  are  indicated  by  dashed  lines.  Using  Figure  1  one  obtains, 
for  Instance,  that  the  undulations  | Aa  j  are  smaller  than  0.01  for  R/ 1^  >  2.0 
in  case  u  =  0.9,  and  for  R/ >  2.5  in  case  a  =  0.5. 

The  relation  (2.3)  also  can  be  expressed  in  terms  of  the  number  N  of 
particles  in  the  averaging  volume.  To  arrive  at  such  an  expression  we  use 
the  approximation 


N  =  C2R/ I^)3 


(2.4) 


(see  Eq.  (3.5)).  Substituting  this  approximation  into  Eq.  (2.3)  one  obtains 


[.An  |  <  2.0  ( 1  -  a)  n2N~2,/3 


(2.5) 


For  spherical  averaging  volumes  Eqs.  (2.3)  and  (2.5)  are  equivalent. 
However,  because  the  latter  also  can  be  used  for  non-spherical  averaging 
volumes,  we  shall  restrict  further  discussions  to  Eq.  (2.5).  Let  | Aoi  | tol 
a  tolerance  for  the  undulations.  Then  one  can  reformulate  Eq.  (2.5)  as  the 
following  condition  for  N 

N  >  [2  (1  -  a  )]3/2  a3  I Aa | _3{2  .  (2.6) 

1  1  tol 


It  N  satisfies  Eq.  (2.6)  then  j  Acx  |  <  ( Aot  |  ^  o  1 ' 

If  i  is  close  to  one,  then  it  is  reasonable  to  replace  the  fixed 
tolerance  by  a  tolerance  level  that  is  a  fraction  of  1  -  a,  for  instance,  by 
the  condition 


Aft  <  k  ( 1  -  a) 


(2.7) 
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with  a  proper  k  <  1.  The  corresponding  condition  on  N  is 


(  i/k) 


'2  ^ 


(2.8) 


One  nuy  combine  both  conditions  on  N  by  requiring  that  N  should  be  larger 
than  the  maximum  of  the  expressions  (2.6)  and  (2.8).  Figure  2  shows  the 
result  for  k  =  0.2  and  various  values  ^f  a  .  According  to  that  figure,  one 
needs  22-110  particles  in  order  to  keep  the  undulations  below  the  tolerance 
levels  ! |  tol  =  0*01  and  k  =  0.2  within  the  range  0.7  S  a  <  1.0.  On  the 
other  hand,  if  one  is  content  with  the  tolerance  level  =  0.03  and  k 

=  0.2.  then  32  particles  suffice  to  represent  reasonable  average  properties. 

One  can  obtain  more  convenient  and  less  sharp  formulas  than  Eqs.  (2.6) 
and  (2.8)  hv  choosing  the  worst  case  conditions,  namely,  a  =2/3  for  Eq . 
(2.61  and  i  =  1.0  for  Eq.  (2.8).  Then  the  combined  condition  on  N  is 


N  >  ma  x 


r  ,2,9/2 


(2/k) 


3/2, 


(2.9) 


The  first  argument  on  the  right  hand  side  of  Eq.  (2.9)  corresponds  to  the 
rightmost  curve  in  Figure  2.  The  undulations  of  averages  of  other  flow 
parameters,  such  as  gas  density,  are  closely  related  to  the  undulations  of 
i.  We  investigate  this  relation  in  Appendix  E  and  show  by  examples  that  the 
amplitude  A;-  of  particle  induced  undulations  of  an  average  gas  property  j>  is 
proportional  to  Ai.  If  the  average  has  a  gradient  and  changes  by  5<j>  along  a 
distance  equal  to  the  diameter  of  the  averaging  sphere, then 

j.'.y|  «  |5?|  |Aa|/(3a)  .  (2.10) 

If  c‘  =  0,  but  there  are  local  particle  induced  inhomogeneities ,  then  the 
undu  i  it  ions  again  are  proportional  to  Aa  and  to  the  amplitude  of  the  local 
inhomogeneities.  In  summary,  undulations  of  a  usually  induce  proportional 
undulations  of  other  flow  parameter  averages. 

Some  consequences  of  these  findings  for  interior  ballistics 
computations  are  discussed  in  Appendix  B,  where  it  is  shown  that  the  mean 
distance  Tm  between  particle  centers  typically  is  as  large  as  three  initial 
propellant  particle  diameters.  According  to  Figure  1,  the  corresponding 
radius  R  of  the  averaging  sphere  should  be  larger  than  six  to  eight  particle 
diameters  (assuming  ) Act  j  1 0  ^  =  0.01).  This  means  that  a  two-phase  theory 

which  is  based  on  averaging  may  not  provide  correct  information  about  radial 
flow  structures  In  Interior  ballistics  because  the  minimum  diameter  of  the 
averaging  sphere  is  approximately  equal  to  the  diameter  of  the  barrel.  Such 
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2.  Minimum  Number  of  Particles  in  an  Averaging  Volume 
for  given  Tolerance  |Aa|tol  an^  ^  “  0*2 


a  theory  still  can  be  used  to  represent  averages  over  cross-sectiona] 
segments  of  the  gun  tube,  that  is,  it  can  be  used  to  calculate  the  interior 
ballistics  core  flow. 

Volume  averaging  smoothes  and  distorts  not  only  the  undesirable 
heterogeneities  ^aased  by  single  particles,  but  also  other  flow  structures, 
particular ly  if  they  have  extensions  less  than  a  diameter  of  the  averaging 
volume.  An  example  of  a  flow  structure  with  a  small  extension  is  the 
boundary  of  a  particle  aggregate.  Let  the  average  gas  volume  fraction  in 

the  aggregate  be  a  and  let  the  aggregate  occupy  the  half-space  z  >  0. 

Conceptually  one  thinks  of  the  boundary  at  z  =  0  as  a  narrow  transition 
between  a  region  with  a  =  1  (gas  only)  and  a  region  with  a  =  a  (gas  particle 
mixture).  The  correctly  defined  a  has,  instead  a  finite  transition  zone 
between  ct  =  3  and  a  =  a  with  a  width  equal  to  a  diameter  of  the  averaging 

volume.  If  the  volume  is  a  sphere  with  the  radius  R,  then  the  transition  is 

approximately  given  by  the  following  function  (see  Appendix  D) 

a  =  1  if  z  <  -R 

ct(z)  =  a  =  1  -  (1^)  (1+|)2  (2-~)/4  if  -R  <  z  <  R  ,  (2.11) 

ot=a  if  R  <  z  . 


This  transition  is  an  idealization,  derived  under  the  assumption  that  the 
undulations  Aa  are  zero  within  the  particle  aggregate,  but  it_  is  not  the 
limit  curve  for  R  >  00 .  (Such  a  limit  is  the  constant  a  *  (l-hx)/2  for  any 
aggregate  occupying  a  half-space).  The  real  transition  curve  undulates 
around  a  basis  that  is  approximated  by  Eq.  (2.11).  The  amplitude  of  the 
undulations  is  bounded  by  Eq.  (2.3)  in  which  u  is  replaced  by  the  basis 
value.  The  approximation  (2.11)  is  quite  good  even  for  moderate  values  of 
R/L,,,,  as  shown  by  the  examples  of  transition  curves  in  Section  4. 

3.  GAS  VOLUME  FRACTION  AS  A  FUNCTION  (F  THE  AVERAGING  VOLUME 

We  consider  the  gas  volume  fraction  a  in  an  averaging  sphere  with  the 
radius  R  and  center  coordinate  vector  £  .  We  assume  that  the  particles  are 
spheres  with  radius  s  and  arranged  in  a  regular  lattice  with  a  lattice 
constant  L.  Let  a  be  the  limit  value  of  the  gas  volume  fraction  for  an 
infinite  averaging  volume.  Then  one  can  define  (see  Appendix  B)  a  mean 
distance  Lm  between  particle  centers  by 


Lm=4s(l-a)“1^.  (3.1) 

The  gas  volume  fraction  depends  on  the  lattice  type,  the  location  X  of  the 
averaging  sphere  with  respect  to  the  lattice,  the  radius  R,  the  radius  s, 


and  on  Che  lattice  constant  L.  Because  a  is  dimensionless,  this  dependence 
can  be  expressed  by  a  function  of  the  following  form 

ot  =  f^S/V  R/ 1^,  s/Luj,  L/  Ljjj)  ,  (3.2) 

whereby  the  function  fj  depends  on  the  lattice  type.  The  ratio  L/l^  Is  a 
constant  for  a  given  lattice  and,  therefore,  it  can  be  included  in  the 
definition  of  the  function  f  ^  .  Also,  because  of  the  definition  (3.1)  the 
argument  s/I^  can  be  replaced  by  a  .  Hence,  one  has  the  following  two 
equivalent  representations  of  the  functional  dependence  of  a  in  a  given 
lattice: 

a  -  f2  (S/I^,  R/Lgj,  s/1^)  (3.3) 

and 


a  =  f3  (fc/Lfc,  R/Lfc,  a)  .  (3.4) 

Eqs .  (3.3)  and  (3.4)  show  that  a  is  a  function  of  five  dimensionless  scalar 
parameters:  the  three  components  of  the  position  vector  Ic/L^,  the  averaging 

radius  R/L^  and  s/L,,,  or  a  .  In  this  section  we  investigate  the  dependence 
of  a  on  the  parameter  R/L^.  We  may  think  of  this  dependence  as  the  result 

of  placing  the  averaging  sphere  at  a  fixed  position  $/Lm  within  a  given 

lattice  and  letting  its  radius  R/L^  vary.  If  R  is  small,  then  one  obtains 

either  a  =  1  or  a  =  0,  depending  whether  the  fixed  center  position  is 

outside  or  inside  a  particle.  As  R/ increases,  a  approaches  the  value 
a  regardless  of  the  center  location.  Figures  3  through  7  show  typical 
curves  of  the  function  a(R/Lm)  for  four  different  lattices.  (The  lattices 
are  defined  in  Appendix  A  and  a  was  calculated  with  the  algorithm  described 
in  Appendix  C) . 

Figure  3  illustrates  the  dependence  of  a  on  R/l^  in  a  square  cylinder 
lattice  with  a  >•  0  .5 .  The  central  graph  shows  two  curves .  One  curve 

corresponds  to  the  center  location  coordinates  ^t/ «■  (0,0,0)  and  starts 
with  a=0,  because  the  point  of  origin  is  a  lattice  point  and  occupied  by  a 
solid  sphere.  The  other  curve  starts  with  a  ■  1  .  It  corresponds  to  the 
center  coordinates  ■  (0.4,  0.2,  0),  which  is  a  point  occupied  by  gas. 

The  upper  graph  in  Figure  3  shows  the  same  two  curves  in  a  larger  scale . 
One  sees  that  the  undulations  of  a(R/l^)  around  the  limit  value  a  have  a 
wave  length  of  about  one  and  an  amplitude  that  seems  to  decrease 

proportionally  to  a  power  of  R/ 1^ .  The  lower  graph  in  Figure  3  shows  the 

number  N  of  particles  that  are  inside  the  averaging  sphere  also  as  a 
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function  of  R/l^.  The  two  curves  corresponding  to  the  two  locations  of  the 
averaging  sphere  are  almost  identical  and  both  approximate  the  cubic  (see 
also  Eq  .  (8.6)  in  Appendix  B) 

N  =  ( ^R/ L^,) 3  .  (3.5) 

Figures  4  through  7  present  the  same  functions  as  Figure  3  but  for 
different  values  of  a  and  for  different  lattices.  One  notices  that  the 
qualitative  as  well  as  quantitative  behavior  of  the  functions  are  very 
similar  in  all  cases. 

Our  goal  is  to  obtain  an  estimate  for  the  amplitude  of  the  undulations 
of  a  about  the  limit  value  a  as  a  function  of  R/ .  To  that  end  one  can 
determine  ror  a  fixed  a  and  R/L^  the  extreme  values  of  a  by  treating 
^/L  as  s  free  parameter  vector.  (One  can  envision  this  as  moving  the 
aveiaglng  sphere  with  a  fixed  radius  in  an  infinite  lattice  until  the  two 
positions  are  found  that  correspond  to  minimum  and  maximum  a.)  The 
calculations  were  done  with  a  simplex  algorithm  and  some  results  are  shown 
in  Figures  8  through  15  .  The  figures  show  that  there  are  differences  in 
detail  between  lattices  and  between  cases  with  different  a,  but  the  general 
trend  of  the  undulation  amplitude  is  rather  uniform  in  all  cases.  This 
trend  is  illustrated  in  Figures  16  and  17.  Figure  16  shows  in  the  left  hand 
graph  the  extreme  negative  deviations  a-a  from  Figure  15  plotted  over  the 
corresponding  values  of  R/I^  in  log, log-scale .  The  right  hand  graph  shows  a 
similar  plot  of  the  extreme  positive  deviations.  Both  graphs  of  Figure  16 
are  combined  into  one  graph  in  Figure  17.  It  is  apparent  from  these  figures 
that  the  negative  and  positive  deviations  behave  similarly  and  that  their 
amplitudes  are  in  general  proportional  to  (R/I^)-^.  Plots  of  the  extreme 
deviations  for  other  lattices  look  very  much  like  Figures  16  and  17 .  A 
combined  plot  of  deviations  in  all  four  lattices  for  a  =  0 .5  is  shown  in 
Figure  18  and  for  a  =  0.9  in  Figure  19.  The  different  symbols  in  these 
figures  indicate  different  lattices  and  different  signs  of  the  extreme 
deviations.  The  dispersion  of  the  plotted  points  shows  that  none  of  the 
considered  lattices  consistently  produces  larger  undulations  than  other 
lattices,  and  that  neither  negative  nor  positive  deviations  are  consistently 
larger  within  the  R/ range  considered.  An  estimated  upper  bound  of  the 
deviations  is  shown  as  a  straight  line  in  the  log , log-plots .  It  has  the 

equat ion 

|*«|  =  0.5  (l^)a2(R/Lm)“2  .  (3.6) 

Eq .  (3.6)  estimates  a  bound  for  the  maximum  amplitude  of  undulations  for  a 
given  size  of  averaging  sphere.  Using  Eq .  (3.5)  one  can  also  express  the 
estimate  in  terras  of  the  number  N  of  particles  within  the  averaging  sphere: 

(Continued  on  page  33) 


;ure  3.  Gas  Volume  Fraction  Dependence  on  Averaging  Radius. 
Square  Cylinder  lattice,  a  *  0.5;  s/l^,  *  0.3969 
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Figure  15.  Extremes  of  Gas  Volume  Fraction  as  Functions  of  Averaging  Radius 
Leap-Frog  Triangular  Lattice,  a  =  0.9;  s/L  =  0.2321 


Figure  17.  Extreme  Deviations  of  Gas  Volume  Fraction 
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Extreme  Deviation  of  Gas  Volume  Fraction  for  All  Lattices 
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If  one  wants  the  undulations  to  he  below  a  tolerance  level  then  Eqs.  (3.6) 
and  (3.7)  can  he  used  to  estimate  the  corre  ponding  minimum  size  of  the 
averaging  sphere  and  minimum  number  of  particles  within  the  sphere, 
respect ivel /.  A  discussion  of  such  estimates  and  corresponding  graphs  are 
presented  i a  Section  2. 


Next  we  discuss  the  validity  of  Eq.  (3.6)  and  (3.7).  They  are  based  on 
sample  calculations  within  the  range  1.0  <  R/Lm  <  4.0  (or  8  <  N  <  512) 

for  a  =  0,5,  0.667  and  0.9,  and  tor  four  different  lattices.  Because  the 
four  lattices  have  quite  different  symmetries  and  maximum  packing  densities, 
one  can  assume  that  the  results  are  representative  for  all  reasonably 

regular  and  uniform  particle  arrangements.  One  can  of  course  construct 
periodic  particle  arrangements  for  which  the  undulations  are  in  excess  of 
the  estimated  limit  at  some  point.  However,  in  order  to  do  so,  one  would 
have  to  stretch  the  concepts  "unifoim"  and  "regular."  The  calculations 
cover  all  a  values  of  practical  importance,  that  is,  all  a  values  over 

0.5.  (The  formulas  (3.6)  and  (3.7)  both  have  the  correct  limit 

A. i  =  0  for  a  =  1.)  We  notice  that  the  minimum  value  of  a  (closet  packing 

of  spheres)  is  between  0.26  and  0.48  for  the  four  lattices  considered.  (See 
Appendix  A.)  Therefore,  values  of  a  below  0.5  are  only  possible  for 
special  arrangements  of  tue  particles,  and  extrapolation  of  the  general 
formulas  to  a  below  0.5,  should  be  done  with  reservations.  The  range  of 
the  parameter  R/ between  1.0  and  4.0  was  also  chosen  by  practical 
considerations.  For  values  of  R/L^  below  one  the  undulations  are  large  and 
Irregular,  indicating  that  the  corresponding  average  is  not  a  good 
representation  of  a  particle  aggregate.  On  the  other  hand,  at  R/Lm  =  4  the 
undulations  have  an  amplitude  of  the  order  10-^,  which  is  sufficient  for 
most  applications.  Extrapolations  of  the  formulas  to  larger  values  of  R/Lm 
likely  will  give  correct  estimates,  but  they  have  not  been  tested  by  sample 
calculations. 


4.  GAS  V OLD  MR  FRACTION  PROFILES 


In  this  section  we  give  an  example  of  gas  volume  fraction  profiles, 
i.e.,  of  the  dependence  of  the  gas  volume  fraction  a  on  the  trajectory  of 
the  averaging  sphere,  and  of  the  smoothing  effect  of  averaging  on  a  flow 
structure.  The  example  is  a  particle  aggregate  with  a  uniform  limit  gas 
volume  fraction  a  that  occupies  the  half-space  z  >  -s  .  (By  a  uniform  limit 
gas  volume  fraction  we  mean  that  the  particles  are  arranged  in  a  regular 
pattern  (lattice)  with  the  limit  value  a  for  that  part  of  any  infinite 
averaging  sphere  that  is  inside  the  aggregate.)  Let  the  particles  be 
spheres  with  radius  s.  We  define  the  boundary  of  the  aggregate  as  the  plane 
z  =  -s,  that  is,  the  tangential  plane  to  the  border  particles  with  a  zero  z- 
coordinate  of  their  centers.  A  naive  representation  of  the  particle 
aggregate  in  terms  of  a  is  by  the  step  function 


if  z  <  -8 
if  z  >  -s 


(4.1) 


However,  this  representation  is  not  correct  because  the  gas  volume  fraction 
is  only  defined  in  terms  of  an  averaging  volume.  (The  gas  volume  is  by 
definition  an  a-fraction  of  an  averaging  volume.)  The  step  function  (4.1) 
does  not  correspond  to  any  finite  or  infinite  averaging  volume  and, 
therefore,  cannot  be  a  representation  of  the  aggregate  boundary.  A  true 
representation  of  the  aggregate  has  to  include  a  size  parameter  of  the 
averaging  volume,  in  our  case  the  averaging  radius  R,  and  one  obtains 
different  transition  curves  between  a  =  1  (outside  the  aggregate)  and 
a  “  a  (inside  the  aggregate)  depending  on  R,  on  the  trajectory  of  the 
averaging  sphere  crossing  the  boundary,  and  on  the  arrangement  (lattice)  of 
the  particles  . 

Figure  20  shows  examples  of  transition  curves  for  two  different 
averaging  spheres  and  two  trajectories  for  each  sphere.  The  spheres  have 
the  radiuses  R/L^  =  1.0  and  R/l^  =  2.0,  respectively,  and  the  two 

trajectories  are  along  the  z-axis  and  along  the  line  x/I^  ■  0.5,  y/I^,  = 
0.5.  The  particles  are  spheres  with  radius  s  and  they  are  arranged  in  a 
leap-frog  triangular  lattice.  (Lattices  are  defined  in  Appendix  A.)  Figure 
20  shows  the  transition  curves  in  two  different  scales.  The  general  slopes 
of  the  transitions  are  about  (1  -a)/(2R/Lm),  but  their  details  depend  on 
the  trajectories.  Inside  the  particle  aggregate  one  observes  periodic 
undulations  with  an  approximate  wavelength  of  1  .5  about  the  limit  value 
a.  The  amplitude  of  the  undulations  is  bounded  by  Eq .  (3.6)  and  it  of 
course  does  not  decrease  as  the  averaging  sphere  is  moved  deeper  into  the 
aggregate . 
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Figure  20.  Gas  Volume  Fraction  at  a  Particle  Aggregate  Boundary 

a  =  0.5;  s/I^  =  0.3969 
R/ L^  =  1.0  and  2.0 
Leap-Frog  Triangular  Lattice 
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Figures  21  and  22  show  transitions  along  the  same  trajectories  and  the 
same  lattice,  but  for  other  particle  sizes,  that  is, for  other  limit  values 
a.  The  general  patterns  of  the  transitions  and  undulations  are  the  same  as 
in  Figure  20,  and  they  are  typical  for  all  four  lattices  considered  in  this 
report  .  The  wavelength  of  the  undulations  was  found  to  vary  between  1  .0 
and  1  .5  L^,  depending  on  the  particle  arrangement  . 

An  approximate  formula  for  the  transition  curve  can  be  derived  under 

the  assumption  that  the  undulations  within  any  subsection  of  the  cloud  are 

zero  .  Then  the  curve  is  (see  Appendix  D) 

1  if  z  +  s  <  -R 

a(z)  =  (  1-0.23  (1-^)  [  l+(  z+s  )/R]  2  [  2-(  z+s  )/R  ]  if  -R  <  z  +  s  <  R  [4.2) 

a  if  R  <  z  +  s  . 


"V 


» 
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\  representation  of  the  particle  agg.  gate  by  the  function  (4.2)  is  better 
than  Kq  .  (  *  .1  )  ,  because  the  former  fakes  into  account  the  averaging 
radius  .  However,  since  undulations  have  been  neglected,  any  real  transition 
curve  will  deviate  from  a(z)  by  undulations  that  are  bounded  by  Eq  .  (3.6) 

where  a  is  to  be  replaced  by  a(z)  .  In  addition  to  the  undulations,  other 

deviations  can  be  present  due  to  a  general  structure  of  the  particle 
arrangement  . 

The  transition  curve  a(z)  with  bounds  computed  by  Eq  .  0.6)  and 

corresponding  extreme  values  of  a  are  shown  in  Figures  23-23  for  a  leap-frog 
triangular  lattice.  The  extreme  values  were  computed  for  fixed  z-values  by 
a  simplex  algorithm  for  which  the  coordinates  x  and  y  were  treated  as  free 
parameters.  Conceptually  this  means  moving  the  averaging  sphere  within  the 
plane  z=const  .  until  the  two  positions  with  extreme  values  of  a  have  been 
found  . 

According  to  Figure  23  the  approximation  a(z)  is  quite  good  in  this 

particular  case.  However,  if  the  particle  size  is  reduced,  as  in  Figure  24, 

then  a  systematic  deviation  becomes  apparent  in  the  upper  part  of  the 
curve.  This  deviation  is  also  visible  in  Figure  23,  that  is,  it  does  not 
disappear  with  increased  averaging  radius.  Such  deviations  were  also 
obtained  for  other  lattices,  and  they  can  be  positive  as  well  as  negative, 
depending  on  the  lattice  and  on  the  orientation  of  the  boundary. 


5.  DISCUSSION  OF  THE  RESULTS 

The  sample  calculations  in  Section  3  and  4  show  that  the  gas  volume 
fraction  a  is  a  function  that  undulates  about  a  limit  value  ci  with  a 
wavelength  1  .3*1^  ( 1^  is  a  mean  distance  between  particle 

(Continued  on  page  42i 
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centers)  and  with  an  amplitude  bounded  by  Eq .  (3.6).  Now  we  discuss 
consequences  of  this  behavior  of  a  for  the  representation  and  computation  of 
a  flow  field.  First,  we  notice  that  in  an  averaged  flow  field  any  flow 
structure  with  extensions  less  than  R  (the  radius  of  the  averaging  sphere) 
is  smoothed  out  and  also  possibly  masked  by  the  undulations.  Second,  the 
minimum  size  of  R  depends  on  the  demands  one  makes  on  the  averages . 
Obviously  one  can  only  consider  an  average  as  a  property  of  a  particle 
aggregate  if  there  are  many  particles  representing  the  particulate  phase 
within  the  averaging  volume.  This  remark  is  quantitated  by  the  estimate  of 
undulation  bounds,  Eq .  (3.6).  As  the  discussion  in  Section  2  shows,  one 
needs  30-150  particles  to  keep  undulations  of  a  below  0.01  and  at  least  30 
particles  if  the  limit  is  0.03.  The  corresponding  minimum  values  of  the 
radius  R  are  2.4*1^  and  1.8*1^,  respectively.  For  the  sake  of  the  following 
discussion  we  assume  that  R  is  chosen  to  be  equal  to  or  larger  than  2*1^. 

The  findings  have  the  follow,  ig  consequences  for  the  representation  of 
a  two-phase  flow  field  by  average  flew  parameters.  Because  one  cannot 
distinguish  between  noise  (particle  induced  undulations)  and  structures 
which  have  extensions  less  than  R,  one  has  to  classify  as  noise  any  flow 
structures  that  extend  less  than  1^.  Consequently,  if  a  flow  field  is 
specified  in  a  mesh,  the  mesh  constant  need  not  be  finer  than  1^ .  Any  mesh 
refinement  can  be  done  by  interpolation  without  loosing  accuracy .  The  same 
applies  to  experimental  data.  If  measurements  are  done  locally  with  a  fine 
mesh  then  the  results  should  be  averaged  over  a  volume  corresponding  to  R  = 
2*1^  at  least.  This  applies,  of  course,  also  to  measurements  near 
boundaries  where  the  averages  again  represent  a  finite  averaging  volume,  and 
not  a  local  point . 

Similar  consequences  can  be  drawn  for  the  computation  of  two-phase  flow 
fields,  for  instance,  by  solving  average  flow  equations.  A  reasonable 
computing  mesh  constant  is  of  the  order  of  L^.  Refinement  of  the  mesh 
constant  below  has  at  best  the  effect  of  interpolation  and,  at  worst,  it 
may  induce  noise  into  the  solution.  In  any  case  such  a  refinement  would 
waste  computing  time.  If  flow  structures  with  extensions  less  than  R  appear 
in  the  solution,  for  instance,  in  mesh  refinement  studies,  then  they  should 
be  considered  as  numerical  artifacts,  because  they  cannot  be  interpreted  as 
average  flow  properties. 

These  considerations  have  an  interesting  result  for  transient  flow 
fields  in  which  a  +  1  in  some  parts  of  the  field.  If  a  approaches  one 
because  the  particle  sizes  are  duced  (for  instance,  by  combustion)  then 
the  condition  R  >  2*  is  not  violated  and  the  particulate  phase  can  be 
represented  by  averages.  If,  however,  a  approaches  one  because  particles 
diffuse  from  the  mixed  phase  region  into  gas-only  regions,  then  the 
condition  R  >  2*  will  be  eventually  violated  and  the  average  model 

invalidated.  A  correct  approach  in  such  a  situation  is  to  model  the 
rarlfied  parts  of  the  particle  aggregate  by  some  other  method  than 
averaging,  for  instance,  by  individually  tracing  the  diffused  particles. 
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In  Appendix  B  we  discuss  the  value  of  the  mean  distance  between 
particle  centers  and  conclude  that  in  typical  interior  ballistics  problems 
Ljh  can  be  as  large  as  three  initial  propellant  particle  diameters. 
According  to  Figure  1  one  then  needs  an  averaging  sphere  radius  that  is 
equal  to  6  to  8  initial  particle  diameters  in  order  to  keep  the  undulations 
below  0.01.  For  some  interior  ballistics  situations  this  means  that  no  flow 
details  in  radial  direction  can  be  represented  by  a  two-phase  theory  that  is 
based  on  averaging,  because  the  minimum  diameter  of  the  averaging  sphere 
approaches  the  diameter  of  the  barrel.  In  such  situations  the  two-phase 
theory  only  can  model  averages  taken  over  a  cross-sectional  segment  of  the 
barrel,  that  is,  the  theory  can  be  used  to  compute  the  interior  ballistics 
core  fLow. 
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APPENDIX  A 


LATTICES 


For  Che  calculation  of  the  gas  volume  fraction  a  one  has  to  specify  the 
size  and  the  spacial  arrangement  of  the  particles.  In  this  report,  the 
arrangement  is  assumed  to  be  in  the  form  of  a  three-dimensional  lattice.  In 
order  to  assess  the  significance  of  the  arrangement,  four  different  lattices 
were  used  for  the  calculations.  These  lattices  are  defined  in  this 
appendix  . 

The  square  cylinder  lattice  is  shown  in  Figure  A.l  .  It  may  be 
constructed  by  arranging  the  particles  in  a  square  pattern  with  a  mesh 
constant  L  in  the  x,y-plane,  and  translating  the  arrangement  in  the  z- 
direction  by  multiples  of  L.  Each  original  square  thereby  generates  a 
square  cylinder.  The  lattice  is  identical  to  the  well-known  cubic  lattice. 

We  identify  each  lattice  point  by  three  integer  parameters,  i,j ,  and 
k .  The  Cartesian  coordinates  of  a  lattice  point  with  the  identification 

(i,.j,k)  are 


x(i  ,  j ,k)  =  iL  , 

y(  i  ,  j  ,k)  =  jL  ,  (A  .1 ) 

z(i,.j  ,k)  =  kL 


The  centers  of  the  particles  are  assumed  to  be  located  at  the  lattice 
points.  Let  s  be  the  radius  of  a  particle.  If  the  square  cylinder  lattice 
occupies  an  infinite  space,  then  the  corresponding  gas  volume  fraction  is 


-4 


(A  .2) 


The  smallest  value  of  a  is  obtained 
possible .  Because  s  cannot  be  larger 
particles,  one  obtains  from  Eq  .  (A. 2) 


if  the  particles  are  as  large  as 
than  L/2  without  intersection  between 


u  -  1  -  a/6  =  0  .476  401  2 
ini  n 


(A  .3 ) 
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The  triangular  cylinder  lattice  is  shown  in  Figure  A. 2.  It  can  be 
constructed  by  arranging  the  particles  in  an  equilateral  triangular  pattern 
with  the  mesh  constant  L  in  the  x,y-plane  and  translating  the  arragement  by 
multiples  of  L  in  the  z-direction.  Each  triangle  thereby  generates  a 
triangular  cylinder.  Again  using  the  integer  triple  (i,j,k)  for  the 
identification,  of  each  lattice  point,  one  obtains  the  following  formulas  for 
the  Cartesian  coordinates  of  the  points: 

x(i,.i,k)  =  iL  +  |  j  (mod  2)  |  i  L  , 

>*0  =  j  L  ,  (A  .4) 

z( i , j ,k)  =  k  L 

In  a  triangular  cylinder  lattice  that  occupies  an  infinite  space  the  gas 
volume  fraction  is 


a 


8  r  (-S',  3 

-  1  T  > 

3/3  L 


(A  .5) 


where  s  is  the  radius  of  the  particles.  The  smallest  possible  value  of  a  is 
obtained  for  the  largest  possible  value  of  s.  Because  s  cannot  be  larger 
than  L/2,  one  obtains 


=  1  -  ti  / (3/3)  =  0.395  400  2 


(A. 6) 


The  leap-frog  square  lattice  is  shown  in  Figure  A  .3  .  It  can  be 
constructed  by  first  arranging  the  particles  in  the  x,y-plane  in  a  square 
mesh  with  the  mesh  constant  L  and  the  square  sides  parallel  to  the 
coordinate  axes.  Then  the  pattern  is  translated  by  multiples  of  L//2  in  the 
z-direction  and  by  multiples  of  L/2  in  the  x-  and  y-directions .  Thus,  each 
square  is  translated  in  a  leap-frog  manner  from  one  z-plane  to  the  next  . 
The  ensuing  lattice  is  also  known  as  "face -centered  cubic  lattice"  with  the 
Lattice  constant  2/T* L.  Identifying  each  lattice  point  by  the  integer 
triple  (1, i,k),  one  obtains  the  following  expressions  for  the  Cartesian 
coordinates  of  the  points: 

x( 1 , j ,k)  =  i  L  +  | k  (mod  2) |  y  L  , 

y ( i J , k  >  =  j  L  +  jk  (mod  2) |  ^  L  ,  (A. 7) 

?-  ( 1  ,  j ,  k )  =  k  —  L 

/T 
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If  a  leap-frog  square  lattice  occupies  an  infinite  space,  then  the  gas 
volume  fraction  Is 


u  =  1  -  (f)3  .  (A  .8) 

3  n  L 

Because  the  maximum  possible  value  of  the  particle  radius  s  is  L/2,  one 
obtains  from  Equation  (A  .8) 


a  ,  =  1  -  it/(3  /2)  -  0  .259  519  5 

rain 


(A  .9) 


The  leap-frog  square  lattice  and  the  leap-frog  triangular  lattice,  which  we 

describe  next,  have  the  highest  possible  particle  concentrations  (the  lowest 

possible  a  ,  )  of  all  three  dimensional  lattices. 
v  min 

The  leap-frog  triangular  lattice  is  shown  in  Figure  A  .4  .  It  can  be 
constructed  by  first  arranging  the  particles  in  the  x,y-plane  in  an 
equilateral  triangular  mesh  with  the  mesh  constant  L,  as  shown  in  Figure 
A .2  .  Then  the  pattern  is  translated  in  the  z-direction  by  multiples  of  /  2/3  L, 
and  in  the  y-direction  alternatively  by  ±(/3V3)L.  Thus,  the  triangular  mesh 
is  shifted  in  a  leap-frog  manner  as  one  proceeds  from  one  z-plane  to  the 
next.  Identifying  each  lattice  point  with  the  integer  triple  (i,j,k),  one 
obtains  the  following  Cartesian  coordinates  of  the  lattice  points: 


x(i , j ,k)  =  1  L  +| j  (mod  2) |  j  L  , 

/3  /? 

y ( i  ,  j , k )  =  j  —  L  +  |k(mod  2)  j  — L  ,  (A  .10) 

z (  i  , .1  ,k)  =  k  /  ^  L 

The  pas  volume  fraction  a  in  an  infinite  leap-frog  triangular  lattice  can  be 
computed  by  the  same  formula  (A  .8)  as  for  the  infinite  leap-frog  square 
lattice.  Consequently ,  also  the  smallest  possible  a  is  given  by  Eq . 
( A  .9 )  .  ,tU  U 

The  formulas  for  a  have  for  all  four  lattices  the  form 

x  =  1  -  A  (s/L)3  ,  (A  .11) 

where  A  is  a  constant  that  depends  on  the  lattice  type.  Because  we  have 
defined  the  lattice  constant  L  In  such  a  manner  that  it  is  the  closest 
distance  between  any  two  centers  of  particles,  the  minimum  value  of  a  is  In 
ill  cases  obtained  by  setting  s  =  L/2  in  Eq .  (A  .11).  Therefore, 


(A. 12) 


where  8max  Is  the  highest  possible  solid  volume  fraction  In  an  infinite 
lattice.  Eliminating  A  between  Eq.  (A.  11)  and  (A.  12)  one  obtains 


a 


1-83  (s/L)3 

max 


(A. 13) 


as  a  general  formula  for  a.  The  value  of  8  depends,  of  course,  on  the 
lattice  type.  Table  A.  1  lists  the  values  of  am^n  and  8  for  the  four 
lattice  types  considered.  The  last  column  of  the  table  contains  the  ratio 
of  the  mean  distance  between  particle  centers  and  the  neighbor  distance. 
That  ratio  is  also  a  characteristic  of  each  lattice.  The  mean  distance  is 
defined  in  Appendix  B. 
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APPENDIX  B 


NEIGHBOR  DISTANCE  AND  MEAN  DISTANCE 

The  density  of  a  particle  aggregate  can  be  characterized  by  various 
parameters:  the  solid  volume  fraction  B,  the  particle  volume  and  the 

number  density,  the  particle  radius  and  a  mean  distance,  the  particle  radius 
and  a  neighbor  distance  in  a  given  lattice,  etc.  In  this  Appendix  we 
discuss  the  concepts  neighbor  distance  and  mean  distance. 

As  neighbor  distance  we  define  the  smallest  distance  between  particle 
centers.  For  the  lattices  defined  in  Appendix  A  the  neighbor  distance  is 
equal  to  the  lattice  constant  L.  However,  the  number  of  neighbors  with  a 
distance  L  is  different  for  different  lattices.  Thus,  in  a  square  cylinder 
lattice  each  particle  has  six  nearest  neighbors,  in  a  triangular  cylinder 
lattice  the  number  is  eight,  and  in  the  leap-frog  lattices  there  are  twelve 
nearest  neighbors  to  each  particle.  One  can  also  construct  regular  lattices 
in  which  the  number  of  nearest  neighbors  is  not  the  same  for  all  particles. 

For  any  of  the  four  lattices  considered  in  this  report  one  can 
calculate  for  given  a  the  corresponding  neighbor  distance  by  solving  Eq. 
(A. 13)  for  L: 

L  =  2s  (B  /(l-a)]1/3=  2s( B  /S)  1/3  .  (B.  1) 

max  '  ^  max  J 


The  values  of  Bmax  are  given  in  Table  A.  1  for  each  lattice.  We  see  from  Eq. 
( B.  1)  that  the  neighbor  distance  is  different  for  different  lattices,  even 
if  s  and  a  are  the  same.  This  is  not  very  convenient  when  computations  are 
compared.  Therefore,  we  also  define  a  mean  distance  1^  that  is  independent 
of  the  arrangement  of  the  particles. 

Let  the  number  of  particles  in  the  particle  aggregate  be  m,  the  volume 
of  each  particle  be  v(s)  =  (4/3)  n  s  ,  and  the  volume  of  the  particle 
aggregate  be  W.  One  may  conceptually  assign  to  each  particle  the  fraction 
W/m  of  the  total  volume  and  represent  this  fraction  as  a  virtual  sphere  with 
the  diameter  L,n.  This  diameter  we  define  as  the  mean  distance  between  the 
centers  of  particles.  A  relation  between  the  solid  volume  fraction  3  and  L^ 
can  be  derived  as  folLows.  By  definition  one  has 


W  =  n  3 
m  A  m 


(B.2) 


and  the  solid  volume  fraction  in  W  is 


6 


(B  .3) 


m  v(s)  _  v(s) 
W  W7m“ 


8  (s/L  )' 
m 


or 


L 

m 


s/8  1/3= 


2s(l-^)-1/3 


(B.4) 


In  Eqs .  (B.3)  and  (B.4)  we  have  used  3  and  a  (the  limit  values)  because  the 
definition  (B.2)  pertains  to  the  whole  particle  aggregate.  We  notice  that 
the  quotient  W/m  is  the  inverse  of  the  number  density  of  particles. 
Therefore,  is  proportional  to  the  cube  root  of  the  inverse  of  the  number 
density.  If  we  consider  a  finite  part  of  the  particle  aggregate,  defined  by 
an  averaging  volume  with  radius  R  and  containing  N  particles,  then  the 
equation  corresponding  to  Eq  .  (B.2)  is  only  approximately  valid: 


|  ( 2R) 3 
N 


I  L3 


(  B  .5 ) 


(The  equation  is  not  exactly  satisfied  because  is  the  mean  distance  for 
the  whole  particle  aggregate.)  Solving  Eq  .  (B.5)  for  N  one  obtains 


N  -  ^~}3  (  B  .6 ) 

m 

as  an  estimate  of  the  number  of  particles  in  an  averaging  sphere.  The 
approximation  by  Eq .  (B.6)  is  quite  accurate,  as  shown  by  the  numerical 

results  in  Section  3.  The  relation  (B.6)  is  of  course  independent  of  the 
particle  arrangement,  i  .e .,  of  the  lattice  type. 

The  advantage  of  the  mean  distance  over  the  neighbor  distance  L  is 
that  the  former  is  independent  of  the  spacial  arrangement  (lattice)  of  the 
particles.  It  is  a  parameter,  like  a,  3,  s  and  the  number  density,  that 
defines  a  general  characteristic  of  the  particle  aggregate.  Because  0  =  1-a 
and  because  Eqs.  (B.2)  and  (B.3)  hold  between  these  five  parameters,  only 
two  of  them  can  be  prescribed  for  a  particular  aggregate.  In  addition  to 
the  two  parameters,  one  may  also  specify  a  structure  (lattice)  of  the 
particm  arrangement,  which  includes  a  definition  of  a  lattice  constant  L. 
The  relation  between  the  lattice  constant  L  and  the  mean  distance  L^,  is 
different  for  different  lattices.  For  the  lattices  defined  in  Appendix  A 
one  obtains  with  smax  =  L/2  from  Eqs.  (B.l)  and  (B.3) 


L/Lm  =  1 


1/3 

max 


( B  .7 ) 
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Hie  value  of  3max  depends  on  the  lattice  type.  Numerical  values  of  the 
ratir  L^/L  are  given  in  Table  A.l,  and  according  to  that  table  the  mean 
distance  in  the  four  regular  three  dimensional  lattices  is  10%  to  24%  larger 
than  the  neighbor  distance. 

Next,  we  estimate  the  mean  particle  distance  for  interior  ballistics 
applications  .  Obviously  any  average  distance  for  the  whole  gas-particle 
system  increases  during  the  firing  cycle,  because  the  available  volume 
increases  due  to  the  motion  of  the  projectile.  Let  a  be  the  initial  gas 
volume  fraction  in  the  gun  chamber  and  sQ  be  the  initial  radius  of  the 
particles.  Then  the  initial  mean  distance  is  according  to  Eq .  (B.4) 

Lmo  =  ZSq^-V'173  *  (B-8) 

Let  zc|lam],er  be  the  length  of  the  chamber  and  travel  be  the  travel  of  the 

projectile  from  its  Initial  position.  The  ratio  of  the  available  volume  WQ 

before  and  W  after  the  projectile  has  moved  by  travel  aPProxlmately 

1  +  z  .  / ?.  .  ,  The  corresponding  relative  increase  of  the  mean 

travel  chamber.  v  6 

distance  is  according  to  Eq .  (B.2)  the  cube  root  of  this  ratio. 

Therefore,  we  have  the  following  general  formula  for  the  mean  particle 

distance  in  a  gun: 


\  =  WW/Wo)1/3 


=  2s  ( l+z  Jz  )1/3(l^t  )  1/3 

o  travel  chamber  o 


(B  .9) 


According  to  this  equation  the  maximum  of  the  mean  distance  is  reached 
during  a  firing  cycle  when  ztravel  is  the  muzzle  value.  Typical  numerical 
values  in  Eq .  (B.9)  are  as  follows.  The  initial  gas  volume  fraction  aQ  is 
between  0.4  and  0.6,  and  z  travel  ^chamber  >  Therefore,  the  maximum  of 

the  ratio  Lm/2sQ  is  for  a  typical  gun  between  2.64  and  3.02.  Differently 
formulated,  the  maximum  of  the  mean  distance  between  particle  centers  is 
for  interior  ballistics  typically  about  three  initial  particle  diameters. 
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APPENDIX  C 


ALGORITHM  FOR  GAS  VOLUME  FRACTION  CALCULATION 

The  gas  volume  In  the  averaging  sphere  with  radius  R  can  be  calculated 
by  first  computing  the  sum  of  the  volumes  of  all  particles  that  are  located 
Inside  the  sphere  and  subtracting  the  sum  from  the  sphere's  volume.  For 
particles  that  are  intersected  by  the  surface  of  the  averaging  sphere  one 
only  has  to  add  to  the  sum  the  volume  of  the  intersection  between  the  two 
spheres.  Next  we  compute  the  intersection  volume.  According  to  Figure  C.  1 
we  have  the  following  relations  if  two  spheres  intersect; 


A  +  a  =  D 


A2  +  h2  =  R2 


(C.  1) 


a2  +  h2  =  s2 


Eliminating  h2  and  solving  for  A  and  a  one  obtains  from  Eq .  (6.1) 


Center  of 
Averagi ng 
Sphere- 


Center  of 
Particle 


Figure  C.  1.  Intersection  of  Spheres 


A  =  yjr  '  D2  +  R2  -  s2) 
a  =  ~  !  D2  -  R2  +  s2l 


( C.  2 ) 


The  volume  of  the  Intersection  Is  the  sum  of  two  spherical  segment 
volumes  .  The  segment  of  the  averaging  sphere  has  the  volume 


vr  ~  j-  tt  (R-A)  [  3h2  +  (R-A)2]  =  ~  tt  (R-A)2  (2R  +  A) 


( C  .3) 


the  volu ae  of  the  particle  segment  is 


Vs  *  tt  (s  -  a)z  (2s  +  a)  , 


( C  .4 ) 


and  the  volume  of  the  intersection  is 


vi  =  Vr  +  vs 


(C.5) 


The  computation  of  the  solid  volume  within  the  averaging  sphere  can  be 
done  as  follows.  First,  using  the  given  coordinates  of  the  center  of  the 

sphere  and  the  specifications  of  the  lattice,  one  defines  a  search  area  in 
terms  of  minimum  and  maximum  values  of  the  integers  (i,j,k)  (see  Appendix 

A)  .  Next,  one  computes  for  each  lattice  point  within  the  search  area  the 
distance  D  between  the  point  and  the  center  of  the  averaging  sphere .  If 
D  >  R  +  s,  then  the  particle  is  outside  the  averaging  sphere.  If  D  < | R  —  s| 
then  on-:  of  the  two  spheres  (averaging  sphere  and  particle)  is  inside  the 

other,  depending  on  which  one  is  bigger.  Then  the  solid  volume  is  equal  to 
the  sraa '  ler  of  the  two  sphere  volumes.  Finally,  if  |R-s|<D<R+s,  then  the 
intersection  volume  is  calculated  by  the  formulas  (C.2)  through  (C.5)  and 
added  t>  the  sura  of  solid  volumes  that  are  located  inside  the  averaging 
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APPENDIX  D 


GAS  VOLUME  FRACTION  AT  A  PARTICLE  AGGREGATE  BOUNDARY 

Let  a  particle  aggregate  with  the  limit  gas  volume  fraction  a  occupy 
the  half-spice  z  i>  0.  Our  goal  is  to  obtain  an  approximate  transition  of 
the  gas  volume  fraction  from  a  =  1  at  large  negative  values  of  z  to  u  “  <s  at 
large  positive  values  of  z.  In  order  to  derive  the  approximation  we  assume 
that  the  gas  volume  fraction  can  be  approximated  by  a  in  that  part  of  the 
averaging  sphere  that  intersects  with  the  half-space  z  >  0.  Under  this 
assumption  u  i  ■;  only  a  function  of  the  coordinate  zQ  of  center  of  the 

averaging  spher  .  CXir  goal  is  to  find  that  function. 

We  notice  that  the  assumption  about  the  approximation  of  a  by  a  cannot 
be  realized  in  any  real  particle  aggregate.  The  undulations  of  a  about  a 
ran  be  reduced  to  arbitrarily  low  levels  if  R/  is  sufficiently  large,  but 

..nly  it  the  whole  averaging  sphere  is  within  the  half-space  z  >  R.  If  the 

averaging  sphere  is  allowed  to  intersect  partly  with  the  half-space  z  >  0, 
then  the  intersecting  volume  cannot  be  assumed  always  large,  even  if_R/*fa  ^ 
l,  and,  theref  >re,the  undulations  are  only  bounded  by  j  Act  j  <  raax(a,  1-a) . 

Let  zc  be  the  z-coordlnate  of  the  center  of  the  averaging  sphere.  The 
volume  of  the  intersection  of  the  sphere  with  the  half-space  z  >  0  is 


0 

.  if 

z  <  -R  , 

c 

*  (R  +  z  '2(2R  -  z  ) 

J  v  C'  c 

,  if 

-R  <  z  <  R 
c 

(D.  1) 

T  r  R 

,  if 

R  <  z 

e 

Because  we  assume  that  a  -  a  in  any  finite  part  of  the  half-space  z  >  0,  the 
gas  volume  fraction  in  the  averaging  sphere  is 


l 


a  +  [  Vr 


V- 

a' 


I  -  (  1  -  t)  V-/V_ 
a  R 


1) 


(D.2) 


where  VR  =  (4/3)  t  R-*  is  the  volume  of  the  averaging  sphere.  Combining  Eqs. 
i  0.  1)  and  (D.2)  one  obtains 


1 


if  z  <  -R 

c 


a  (z  ) 

c 


l-(l-a)  U  +  z  /R)2  (2  -  z  /R )/ 
c  c 

a 


if  -R  <  zc<  -R  ,  (D.3) 

if  R  <  zc 


The  limit  of  the  transition  curve  for  R/l^*  “  is  simply  a=(l+a)/2. 
This  is  a  nice  Illustration  to  the  smoothing  effect  of  averaging:  If  the 

averaging  volume  is  made  so  large  that  the  particle -induced  heterogeneities 
are  zero,  then  also  all  flow  structures  (aggregate  boundary  in  the  present 
example)  are  reduced  to  constants. 
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APPENDIX  E 


UNDULATIONS  OF  AVERAGE  FUNCTIONS 

Let  |  be  a  local  gas  property,  for  Instance,  density 
corresponding  average  <J>  is  defined  by 


gas  V 

gas 


Then  the 


(E.l) 


where  VT.  .  is  that  part  of  the  averaging  volume  V  which  is  occupied  by 


gO,s  -  ■ 

gas.  if  v  is  constant,  then  obviously  =  41  for  any  positive  value  of  V  g  = 
a  V  .  [f  £  is  not  constant,  then  the  particles  in  V  influence  the  value  of 

the  average  ■£  ,  that  is,  depends  in  such  a  case  on  a.  Particularly,  if 

a  undulates,  then  also  <j>  undulates.  We  estimate  the  amplitude  A<f>  of  the 
undulations  for  a  local  function  <j>  with  a  constant  gradient.  We  assume 
without  loss  of  generality  that  only  the  x-coroponent  of  the  gradient  is  non¬ 
zero  and  express  <{>  by 


4,(x,y,z)  =  +  (x-xq)  $x 


(E.2) 


with  constant  $  ,  xQ  and  <t>  ^ .  The  corresponding  average  value  is  obtained  by 
substituting  the  expression  (E.2)  into  the  integral  (E.l).  Let  the  particle 
volumes  v  be  all  equal  and  small  compared  to  the  averaging  volume  V,  let  x 
denote  the  x-coordtnate  of  the  center  of  V,  and  x^ ,  i=l,...,N,  denote  the 
centers  of  the  N  particles  that  are  inside  the  averaging  volume  V.  Then  the 
average  value  of  $  is  approximately  equal  to 


ar--=-  V  q> 
V-L V  1  T  o 


+(x-xo)^x]  -  V  V  [*o+  (xr  X0HJ} 


V-Nvji^  -(x-x  H  ]  -  v  J  l  (x  -x)l 
' '  o  o  x1  x  ~  i  ' 


(E.3) 


+■  (x-x  )  <{>  -  K  (l-<i)/a. 

o  o  x  x 


where  r,  is  the  average  deviation  of  the  particle  positions  from  the  center 
of  V  . 

Next  we  relocate  the  averaging  volume  to  the  position  x+Ax.  The 
relocation  affects  the  average  value  of  |  in  two  ways.  First,  the  local 
function  is  averaged  over  a  different  region  and  the  positions  of  the  N 
particles  have  changed  relative  a  V.  According  to  Eq  .  (E.3)  this  changes 


71 


the  average  by 


Ax*$  (  1  +  (  l-a)/a) 
x 


(E.4) 


A  second  effect  is  due  to  a  possible  change  of  the  number  of  particles 
within  V.  Lst  AM  be  the  number  of  particles  that  have  been  added  by  the 
relocation  and  Am  be  the  number  of  particles  that  have  been  lost  by  the 
relocation.  The  total  number  of  particles  within  V  then  chafes  by 


=  AM  -  \m 


(E.5) 


and  the  gas  volume  fraction  changes  by 


-■  AM  =  -  (  1-a )AN/N 

V 


(E.6) 


The  average  value  of  <$>  at  x+Ax  is,  with  the  same  approximation  as  Eq .  (E.3), 


<t>  ( x+A x )  =  ■  ,  a \ +  (x+Ax-x  )£  }  - 

V-(N+AN)v  1  L  o  o  xJ 

P 

N  AM  Am 

-  v  [j;  + 1  - 1  )  ;;o+  (xi  -  xo> 


V--(^AN)v-  { [  V-(  N+A  N)  v]  +  (x+Ax-x^*  J 


(E.7) 


N  AM  Am 

~  +  1  ~  l  )  [xi  "  x  ~  Ax] 


>  +  (x+Ax-x  )  4> 

o  ox 


AM  Am 

£x[?-Ax  +  ^  [I  -  l  )  (x1  -  x  -  Ax)] 


We  consider  the  last  terra  in  Eq.  (E.7)  and  observe  that  the  particles  are 
added  and  removed  at  opposite  sides  of  the  averaging  volume.  If  V  is  a 
sphere  with  radius  R, then  one  may  estimate  the  average  x-coordinates  of  the 
added  particles  by  x^«  x  +  2R/3  +  Ax/2  and  the  corresponding  average  of  the 
coordinates  of  the  lost  particles  by  x^°  x  -  2R/3  +  Ax/2  .  Then  the  last 
term  in  Eq.  (E.7)  is 


j  AM  j  fim 

-  I  (x^x-Ax)  -  ^  l  (x1-x-Ax) 


4  AM(2R/3  -  Ax/2)  +  4  Am(2R/3  +  Ax/2)  ^ 
N  N 


(E.8) 


=  4  4  R(A  M  +  Am)  -  4  (A  M  -  Am) 

N  3  N  2 

Next  we  express  AM  and  Am  in  terms  of  a  and  Aot .  The  region  of  space  newly 
covered  by  the  relocation  of  V  has  the  volume 


...  „2  3Ax 

AV  =ttRAx  =  V  — r— 
4R 


(E.9) 


and  an  equal  volume  is  abandoned.  Let  a  ^  be  the  average  value  of  a  in  the 
added  region  and  be  the  average  value  in  the  abandoned  region.  Then  the 

numbers  of  added  and  lost  particles  are 


A  M  =  ( 1  -a  )  — 
M  v 


V  3Ax  M  3Ax 


(E  .10) 


^  Av 
(l-o  )  — 

m  v 


^ra  „  3Ax 

~nr N  ~ar 


Substituting  these  values  into  Eq  .  (E.8),  one  obtains  the  expression 


1  -a  +1  -a  a  ,,-a 

M  m  .  ,  3A  x  Mm 

- - - - A  V  +  - - 

2(l-t)  8R  L-a 


(E.ll) 


It  the  gas  volume  fraction  is  uniform  (a «  a  »  a),  then  this  expression  is 

M  in 

>'qual  to  Ax.  We  assume  that  a  is  not  uniform,  particularly  that  a  is 
increased  by  Aa  due  to  an  excess  of  added  particles  over  the  average  A  M  in 
Kq  .  (K.10).  We  express  this  assumption  by  modifying  Eq .  (E.10)  to 


A  M  =  (l-a)  —  +  AN 
v 


Am  =  ( 1  -<a  ) 


N  —  Ax  +  AN 
4R 


N  4RAX 


(  E  .1 2  ) 
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v  v  r’*  -  *r 


& 


I; 


< 


► 


Then  the  expression  Eq .  E.8)  becomes  equal  to 


Ax  +4PaN  -Iax  {AN 


(E  .13) 


N  3 


N 


Replacing  the  last  term  in  Eq .  (E.7)  by  this  expression  we  obtain 


(x+Ax)  -  £o+  (x+Ax  -xq)4>  x  - 


1-q  ~  f  ,  A  N  /-2R  .  Axu 


or,  using  Eq .  (E .6) 


$(x4Ax)  =  +  (x  -Wx-xn) 


O  x 


1-a  r  r f  _  Aa  2R  f  .  _  3 
a+Aa  x  1-a  3  2 


3  —  )] 
2R;J 


The  difference  between  Eqs .  (E.3)  and  (E.15)  is 


i  (x+Ax)  -  <t>  (x)  =  A  x  p 


(E  .14) 


(E.15) 


(E  .16) 


1 


K  Aa  [ 1  “ 


a+Aa  yx  3 


3  Ax  1-a 

2  2R  a  2RJ 


Hence  the  change  of  the  average  function  $  due  to  a  change  of  a  is 


"1 


y  Aa  r  ,  j  ux  i  -v.  ,  <,  ■ 

C  =  "2R<tx  3 (a  4Aa~)  f  1  ~  2  2R  a  3  2R- 


3  Ax  .  1-a  ,  ? 


(E  .17) 


In  order  to  estimate  the  amplitude  of  the  undulations  of  <p  we  set  Ax  equal 
to  about  one  fourth  of  the  wave  length  of  the  undulations  of  a  •  According 
to  Section  4  the  wave  length  is  between  1.0  and  1.5  1^,-  By  setting  Ax 
1^/3  in  Eq .  < E  .17)  we  obtain 
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1 


(E.  18) 


-  2  R  41 


Aa 


x  3(a+Aa) 


r  1  -  _™  +  ±2. 3  M 

L  4R  ct  2RJ 


The  quantity  1 2R4  |  is  the  maximum  change  of  the  local  gas  property  4  along 
a  diagonal  of  the  averaging  volume.  According  to  Eq.  (E.  16)  this  is  also 
the  change  of  the  average  function  4)  if  undulations  are  not  present.  let 
that  change  be  denoted  by  641.  The  term  L^MR  in  Eq.  (E.  18)  can  be  assumed 
less  than  0.25  or  smaller.  Also,  |s|  «  2R,  so  that  the  two  terms  involving 
1^  and  i  can  be  neglected  in  our  estimate  of  C.  If  we  also  assume 
| Act )  <<  ot,  then  Eq.  (E.  18)  yields  the  following  estimate  for  the  amplitude  A4 
of  the  undulations  of  <j>: 

A*  -  |<5<t>  |  (Aa  j /(3ot)  .  (E.  19) 


A  corresponding  estimate  with  I-a  replacing  a  can  be  obtained  for  the 
average  properties  of  particles,  which  are  defined  by 


1 


r 


v 


4  dV 


parte.  V 


parte. 


(E.20) 


where  Vpartc  is  that  part  of  the  averaging  volume  which  Is  occupied  by 
particles  and  4  is  a  local  particle  property. 


Undulations  of  a  also  can  affect  such  average  functions  which  have  a 
zero  gradient  (iS<|>=0  ),  if  the  local  property  J  is  affected  by  the  presence 
of  particles.  One  example  of  such  a  situation  Is  a  particle  aggregate  that 
noves  with  velocity  u  through  a  quiescent  viscous  gas.  Then  the  average  gas 
velocity  in  ^a  large  averaging  volume  is  not  zero  but  approximately  a 
fraction  of  u.  The  size  of  the  fraction  depends  on  the  thickness  of  the 
boundary  layer  around  each  particle  and  on  the  number  of  particles,  that  is, 
on  ,x .  A  similar  example  is  the  average  temperature  of  a  gas  In  which  one 
places  particles  with  a  different  temperature.  We  notice  that  in  the  latter 
example  the  average  gas  temperature  undulates  with  an  amplitude  proportional 
to  the  undulations  of  a,  whereas  the  average  particle  temperature  is 
affected  only  by  higher  order  terms.  We  demonstrate  this  remark  by 
considering  the  following  simple  model  in  which  the  gas  has  a  constant 
property  4  j,  in  a  boundary  region  around  each  particle  and  a  constant 
property  4Q  elsewhere  in  the  flow  field.  The  volume  of  the  boundary  region 
around  a  particle  we  denote  by  ve .  (For  a  spherical  particle  with  radius  s 
and  boundary  region  thickness  As  «  s  one  has  e  =  3As/s.)  Then  the  average 
gas  property  4  is 
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<»  =  V=N^~  ^V"NV  -  Nv  e)  4?o  +  Nv  e^] 

■  —  [(a  -  (1-a)  £ )  <J>  +  (1-a)  g<(>  .  ]  ,  (E.21) 

Ot  L  O  D 

=  +  (£.  -  <J>  )  e  ( 1-a )/a 

o  bo 

Recause  the  average  <{>  depends  on  a,  undulations  of  a  will  cause  undulations 
of  $.  For  small  undulations  Aa  <<  a  one  obtains  from  Eq.  (E.21)  the 
estimate 


A* 


)  e  A a/a 


2 


(E.22) 


that  is.  A#  Is  proportional  to  Aa.  The  average  particle  properties  are  not 
affected  by  Aa  under  similar  conditions.  Let  $  be  a  constant  particle 
property  In  a  boundary  region  within  each  particle,  and  be  a  constant 
particle  property  elsewhere  in  the  particles.  The  average  ip  then  is 

*  =  n7  UNv  “  Nv  c)  *0  +  NV  E 

=  [  (  !-«  -  (l-o)e)  +  ( l“t»)e<f fc]  (E.23) 

"  *o  +  (*b  '  *o)e 


Because  the  average  i|>  is  independent  of  a ,  it  is  not  affected  by  undulations 
of  a . 


To  summarize,  we  observe  that  if  particle  induced  undulations  of 
average  flow  properties  do  occur ,  then  they  likely  are  proportional  to  the 
undulations  Aa  of  the  gas  volume  fraction. 
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